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Abstract 

We consider the problem of metastability for a stochastic dynamics with a parallel updating 
rule with single spin rates equal to those of the heat bath for the Ising nearest neighbors 
interaction. We study the exit from the metastablc phase, we describe the typical exit path and 
evaluate the exit time. We prove that the phenomenology of metastability is different from the 
one observed in the case of the serial implementation of the heat bath dynamics. In particular 
we prove that an intermediate chessboard phase appears during the excursion from the minus 
metastable phase toward the plus stable phase. 



*AMS 1991 subject classification: 60K35; 82B43; 82C43; 82C80. Keywords and phrases: stochastic dynamics, 
probabilistic cellular automata, metastability, low temperature dynamics. 



1 



1. Introduction 

Metastable states arise when a physical system is close to a first order phase transition. If the 
system is prepared in the metastable phase, it takes an extremely long time to reach equilibrium. 
In physical experiments it is seen that if the system is not suitably perturbed it remains forever in 
the metastable phase |PLj . 

A rigorous description of metastability cannot be formulated in terms of the standard equi- 
librium statistical mechanics: dynamical models must be considered [CGOVj . The case of the 
stochastic serial dynamics has been discussed, for instance, in |KOl INST| IS]: at each step of time 
one of the spins on the lattice is updated with rates satisfying the detailed balance condition. In 
this set-up it has been seen that starting from the wrong metastable phase, the time needed by 
the system to exit the metastable state, namely the exit time, is exponentially long in the inverse 
of the temperature. Moreover, the exit time is exactly the time needed to see a sufficiently large 
droplet, namely the protocritical droplet, of the stable phase in the metastable background. Hence, 
the equilibrium is achieved via the nucleation of such a protocritical seed. 

It is rather natural to ask oneself in which sense these results depend on the dynamics. In 
this paper we consider a dynamics in which simultaneous spin flips |BCLS| IC] are allowed: the 
single spin flipping rates are those corresponding to the two dimensional nearest neighbors Ising 
interaction. More precisely we study the metastable behavior of a Probabilistic Cellular Automaton 
[Rl ISt] which is reversible with respect to a Gibbs measure derived by an Hamiltonian with four 
body interaction. We show that the exit path from the metastable phase to the equilibrium changes 
dramatically, with respect to the serial implementation of the heat bath dynamics, in particular 
the system visits an intermediate metastable phase before reaching the equilibrium. This is not 
surprising, indeed, as it will be pointed out throughout the paper, there exist many deep differences 
between the evolution of the system under a serial and a parallel dynamics. 

We focus, now, on what we consider the most relevant novelty appearing in the study of metasta- 
bility for parallel dynamics: in Glauber dynamics the system can jump between configurations 
differing at most for one spin, such pairs of configurations are called neighboring configurations. A 
connected domain is a subset of the configuration space such that for any pair of states it is possible 
to find a sequence of pairwise neighboring configurations of the domain joining the two states; the 
system, during its evolution, can visit the whole connected domain without exiting from the domain 
itself. In order to exit a connected domain, the system must necessarily cross its external boundary, 
that is the set of configurations not belonging to the domain, but having a nearest neighbor inside 
it. This sort of "continuity" property is the key property in estimating the exit time, that is in 
establishing the minmax between the metastable and the stable states, namely the minimal energy 
barrier bypassed by any path joining the metastable to the stable state. 

Continuity is absent in the case of PCA's: any configuration is connected to any other, a path 
joining the metastable to the stable state is an arbitrary sequence of configurations starting with 
minus one and ending with plus one. The lack of continuity forces us to develop techniques to 
estimate the energy cost of any direct jump from a subcritical to a supercritical configuration. 

The paper is organized as follows: in Section [2] we define the model. In Section [3] we state 
our results: we first characterize the stable configurations (fixed points for the zero temperature 
dynamics, that is the typical droplets of the right phase plunged into the sea of the wrong phase); 
then we study the tendency to grow or to shrink of such droplets; finally, we construct the subset of 
the configuration space visited by the system in the metastable phase (description of the fluctuation 
around the metastable state) and, via a detailed description of the escape path, we estimate the 
exit time. In Sections [4] and [5] we, finally, prove the Theorems and the Propositions. 
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2. Definition of the model 

In this section we define our model, namely a Probabilistic Cellular Automaton reversible with 
respect to a four body hamiltonian. 

2. 1 . Preliminary definitions 

Let A be a two-dimensional torus containing I? lattice sites, i.e., A C Z 2 is a square containing 1? 
points and having periodic boundary conditions. Let d : (x, y) G A x A — > d(x, y) G [0, +oo) be the 
euclidean distance on the lattice A. For any X, Y C A we define d(X, Y) := mf x <=x,yeY d(x, y)- 

We say that x, y G A are nearest neighbors iff d(x, y) = 1. We say that the set X C A is a cluster 
iff for any x, y G X there exist x±, . . . , X}. G X such that x\ = x, x^ = y and for any i = 1, . . . , k — 1 
the two sites Xj and Xj + i are nearest neighbors. 

Given two integer numbers m > I > 1 and x G A we denote by R x ,£, m a rectangle on the dual 
lattice A + (1/2, 1/2) with side lengths t and m and such that x is the first site of A inside the 
rectangle in lexicographic order. We denote by Ri,m := {x G A : x is inside Re, m } the interior of 
Rem- We will drop x from the notation when it will be not necessary to specify the location of 
the rectangle on the lattice. We say that two rectangles R x ,e,m and R x ',e',m' are interacting (resp. 
non-interacting) iff d(R x ^ )Tn , R x '/' iTn ') < 2 (resp. > V5). 

We associate a spin variable o~(x) = ±1 to each site x G A; the space {1, — 1} A of configurations 
is denoted by S. If a G S and X C A we denote by ax a configuration such that ax{x) = o-{x) for 
any x £ X and o~x{x) is arbitrary for any x G A \ X. 

2. 2 . Definition of the dynamics 

Let o~(x) = ±1, for any x G A, be a spin variable and let 

H^ h (a) ■== H^(a) := - £ (r(x)cr(y) - J>(*) (2.1) 

be the Ising nearest neighbors interaction, with the first sum performed over all the nearest neighbor 
pairs, a G S and h G K. 

Let us introduce the discrete time variable n = 0, 1, . . . and denote by a n the system configu- 
ration at time n. All the spins are updated simultaneously and independently at every unit time; 
the conditional probability that the spin at site x takes value a G {— 1,+1} at time n, given the 
configuration at time n — 1, is 



cxp {-f3H {I) (a, K_i) MM )} 

p x {a\a n -i) :- 



exp {-^iT( J )(a, (o- n _i) A \ {z} )} + exp {-/3#( 7 )(-a, (cr n _i) A \ {a; })} 
1 1 



(2.2) 



!, / or l q - [1 + atanh/? (*) + &)] 

1 + exp{-2/3a(5 (Tn _ 1 (x) + /i)} 

where ±a, (cr n _i) A \{ x } are the configurations equal to cr n „i on A \ {x} and to ±a on {x}, 

Sr(a:) := J] a ( y ) 

yGA: d(x,j/)=l 

for any a G 5 and x G A. The normalization condition ^(ojcn-i) + p x (— a\a n -\) = 1 is trivially 
satisfied. Thus the time evolution is defined as a Markov chain on S with non-zero transition 
probabilities P\(rj\a) given by 

P A (ri\o-)=P A (o-,ri):=l[ Px (ri(x)\a) V(7,r/G5 . (2.3) 
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It is straightforward [Dj that the above Probabilistic Cellular Automaton is reversible with respect 
to the Gibbs measure v\(cr) := exp{— H\(a)} /Z\ with Z\ := X^es ex P{ — Ha(v)} an ^ 

H P / l (a) = H(a) := -0h a(x) - log cosh [0 (S a (x) + h)} . (2.4) 

In other words the detailed balance condition 

P A (a, ??) exp{-H A (a)} = P A (r/, a) exp{-H A ( V )} (2.5) 

is satisfied for any a,ij £ S. The interaction is short range and it is possible to extract the potentials: 
for any a 6 5 we can write 

H(a) - const = —J ^ a(x) - J (<)) ^ <r(x)<T(y) - J <(<)>> ^ o"(x)o-(y) 

x£A ((xy)) ({(xy))) 

(2.6) 

-J^o-(x)cr(y)cr(z) - J ^ a(x)a(y)a(w)a(z) 

Xyz Qxywz 

where the five sums are respectively performed over all the sites in A, the pairs of next to the nearest 
neighbors, the pairs of sites at distance 2, the three site clusters composed of two consecutive not 
parallel pairs of next to the nearest neighbor sites and, finally, over the four site diamond shaped 
clusters. The even coupling constants are 



l 1 _cosh/?(4 + /i)cosh/?(4-/i) p-^oo a 1 
cosh 2 (/3/i) 



J «» = 2 J «<»> = s lo S 7T77^~\ ~ P~-Ph 



j 1 cosh 0(4 - h) cosh 6 (0h) cosh 0(4 + h) p->oo 1 3 
~16° S cosh 4 0(2 + h) cosh 4 0(2 - h) ~ r + 8 

while the odd ones are 

1, cosh 2 0(2 + h) cosh 0(4 + h) p^oc 5 

J = 0h + - log tt— — ~ -0h 

4 & cosh 2 /?(2 - h) cosh /3(4 - /i) 2 K 



(2.7) 



1 , cosh 2 0(2 -h) cosh /3(4 + /t) /3^oo 1 a . 

J = — log ~ — 0h 

16 & cosh 2 0(2 + h) cosh 0(4 - h) 8 



(2.8 



T/ie energy and the zero temperature phase diagram 
The definition of ground states is not completely trivial in our model, indeed the hamiltonian 
H A depends on 0. The ground states are those configurations on which the Gibbs measure v A is 
concentrated when the limit — > oo is considered, so they can be defined as the minima of the 
energy 

E\(a) E(a) := hm = a(x) - £ |S CT (x) + h\ (2.9) 

uniformly in a € 5. Notice that it is possible to write H A (o~) = 0E A (o~) + o(exp{— 0c}) for some 
positive constant c depending on a. 

We consider, now, the case h = 0: E A (a) = — YlxeA I^O^)!- ^ is rather clear that there exist 
four coexisting minima +1, — 1, C e ,C° € S: 

+ l(x) = +l, -l(x) = -l, C e (x) = (-l) Xl+X2 and C°(x) = (-If 1+X2+1 (2.10) 
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Fig. 2.1: The zero temperature phase diagram in the plane h a -h e . The four states +1, — 1, 

C° and C e coexist on the solid line whose ending points are (—4,4) and (4,-4). Each 
dotted line is the boundary between two regions with different ground states coexisting on 

the line itself. 



for all x = (xi,X2) £ A. Notice that C e and C° are the chessboard configurations with plus spins 
respectively on the even and odd sublattices. We define C := {C°, C e }. 

Now, we wonder what happens when h ^ 0: a full description of the zero-temperature phase dia- 
gram requires the introduction of a staggered magnetic field. We consider the new zero-temperature 
energy 

E^' h "(a):=-Y t hMx)-Y / \Sa^) + hx\ , (2.11) 

xeA xeA 

where h Q , h e £ R and h x = h Q (resp. h x = h e ) if x belongs to the odd (resp. to the even) sublattice. 
A simple calculation gives the energy of the four zero-field ground states: 

E ho,h e{+l) = _[Aj [K + ^ + | 4 + ^| + | 4 + h ^ and E ho,h e ^ = El h °'- h *{+l) 

. (2.12) 

E*°' h *(C ) = -^ [ho-h e + \4-h \ + \4 + he\] and E^ he {C e ) = El h °^ he {C°) 

By comparing the four expressions (|2.12p one obtains the zero-temperature phase diagram in Fig. 
12.11 We note that on the line h Q = h e = h, depending on the sign of the magnetic field the ground 
state is either +1 or — 1; but at h = there are four different coexisting ground states. 

2.4- Heuristic description of the low temperature phase diagram 

In this section we give a heuristic argument showing that at finite, but very low, temperature the 
structure of the phase diagram is not changed. More precisely the argument suggests that at h = 
the four states +1, -1, C° and C e still coexist [KVllV]. 

At finite temperature ground states are perturbed because small droplets of different phases 
show up. The idea is to calculate the energetic cost of a perturbation of one of the four coexisting 
states via the formation of a square droplet of a different phase. If it results that one of the four 
ground states is more easily perturbed, then we will conclude that this is the equilibrium phase at 
finite temperature. 

A simple calculation shows that the energy cost of a square droplet of side length n of one of the 
two homogeneous ground states plunged in one of the two chessboards (or vice versa) is equal to 
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8n. On the other hand if an homogeneous phase is perturbed as above by the other homogeneous 
phases, or one of the two chessboards is perturbed by the other one, then the energy cost is 16n. 

Hence, from the energetical point of view the most convenient excitations are those in which 
a homogeneous phase is perturbed by a chessboard or vice versa. Moreover, for each state 
— 1, +1, C e , C° there exist two possible energetically convenient excitations: there is no entropic 
reason to prefer one of the four ground states to the others when a finite low temperature is con- 
sidered. This remark strongly suggests that at small finite temperature the four ground states still 
coexist. 



3. Results and heuristics 

We pose, now, the question of metastability: let h be positive and small; we prepare the system 
in the starting configuration gq = — 1 and we try to estimate the first time at which the system 
reaches +1. 

The two chessboard phases coexist at h = with the minus and the plus phase: it is natural to 
wonder if these phases play a role during the escape from the minus metastable phase toward the 
plus stable phase when the external magnetic field is positive and small. 

The main feature of PC A models is that the system can jump from any configuration to any 
other, in contrast with what happens in serial Glauber dynamics, where transitions are allowed 
only between configurations differing at most for one spin. We remark that in this model the single 
spin flip is not a local event, in the sense that its probability depends on all the spin of the lattice. 
Indeed, given X C A we denote by a x the configuration obtained by flipping in a all the spins at 
sites x G X; if X = {x} for some x G A, then by abuse of notation we will denote a x = = a x . 
Now, by (|2.3p we have that 



PA(cT,a x )=l[p x (a*(x)\a) ]J p y {a{y)\a)=\{p x {-a{x)\a) ]J Py (a(y)\a) , (3.1) 
x&x yeA\x xex yeA\x 

that is the probability to flip the spins inside X depends also on the probability that spins outside 
X are not flipped. Notice that this is true even if \X\ = 1, namely if only one spin is flipped. 

3. 1 . Stable configurations 

First of all we characterize the stable configurations of the system, namely those configurations 
(j£5 such that Pa(o~, a) — > 1 in the limit (3 — > oo. Equivalently, a G S is a stable configuration if 
and only if Pa(cj, n) — > in the limit (3 — > oo for all 77 G S \ {a}. 

We discuss, now, the possible single spin events. In Table [T] we consider a site x and we draw 
all the possible configurations in a five spin cross centered at x. The probability p x (+l\a) to see 
+1 at x is evaluated (we recall p x (—l\o~) = 1 — p x (+l\o~)). From TableHJit is clear that in the limit 
(3 — > 00 the probability associated to a single spin event is either one or zero, in the sequel we will 
respectively say high and low probability events. By (|3.ip it follows that the same limiting behavior 
is valid in general for any transition P\(a,o~ x ) with a G S and X C A; in this sense PCA's are a 
generalization of deterministic Cellular Automata. We remark that for any a G S there exists a 
unique configuration n G S such that the transition a — > ij happens with high probability, that is 

-Pa(C) 7 ?) " — -* 1- We note, moreover, that i] = Ta, where T : a G S — ► Tcr G S is the map such 
that for each x G A 

To-(x) := 



if p x (a x (x)\a) P - 



a(x) if p x (a x (x)\a) P - 



(3.2) 
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1 ~ „-2/3(4-/i) _ , 1 p -2f3(2-h) 

l +e 2/3(4-h) — e X -T 1+e 20(2-h) — e 



+ + 



+ 

+ x + I ~ 1 _ e -2/3(4+h) 

' ^ l +e -2/3(4+h) — 1 e 

+ 



Tab. 1: Probabilities for single spin events: probability to sec +1 at site x at time t, with 
the neighboring configuration at time t — 1 drawn in the picture. 



that is at each site we do the right thing in the sense of following the drift. We can say that <r£5 
is a stable configuration iff cr = Tcr. 

In order to characterize the stable states of the model we need few more definitions: let C £ C = 
{C°, C e }, we denote by Sc C S the set of configurations with a well defined sea of chessboard C. 
Similarly we define <S+i C S and we set Sc ■= Sc°USc e ■ More precisely, for each a G {— 1, +1}, 
for each a £ S a there esists a percolating cluster X C A such that crx = qi and ax = (T n a)x 
for all n > 1; for each a G {C e , C" 5 }, for each a £ S a there esists a percolating cluster X C A such 
that ax = ax arid ax = (T 2n <7)x for all n > 1. 

Proposition 3.1 A configuration a G 5_i is stable for the PCA \2. 3\) iff a(x) = +1 /or a// t/ie 
siies x inside a collection of pairwise non-interacting rectangles of minimal side length i > 2 and 
a(x) = — 1 elsewhere. A configuration a G <S+i is stable iff a = +1. There is no stable configuration 
a £ S c . 

In other words we can say that the only not trivial stable states are configurations with well 
separated rectangular droplets of pluses inside the sea of minuses. The Proposition 13.11 follows from 
[rTST] and Lemma E2J 

Lemma 3.2 A configuration a £ S is stable for the PCA \2. 3\) iff 

p x {a x {x)\a) ^ Vx G A 

Proof of Lemma \3.2\ Suppose p x (a x (x)\a) - — ? for all x £ A: let n £ S\{a}, there exists X C A 
and X/fl such that rj = a x ; thus, by equation (|3.ip one has 

PaM) = Pa(<7 i( 7 x )= n^(^Wk) II Py(°(y)W) ^ o 

xex ye\\x 

Suppose a is a stable configuration: P(a,a) — ► 1 in the limit — > oo, (|2.3p and the normalization 
condition pa;((T(x)|cr) +p a ;(cr a: (a;)|o") = 1 imply the statement. □ 

3.2. Stable pairs and traps 
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The configurations in which our system can be trapped are not exhausted by the stable configura- 
tions. Indeed, let a G S and 77 = Ta / a the unique state reached with high probability starting 
from a. If it were T77 = a, then the system would jump back and forth from a to r\ with probability 
going to one in the zero temperature limit; the system would be trapped into a two state loop. 
Given a,rj G S and a 7^ 77, we say that they form a "stable pair" iff 77 = Ta and T77 = a. The two 
chessboard configurations C° and C e are a simple example of a stable pair. 

We discuss two important properties of the stable pairs. From the detailed balance condi- 
tion (j2.5p it follows that if cr, 77 6 S form a stable pair, then they have the same energy, namely 
E A (a) = Ea.{tj). Indeed, from (j23j) and (JZJ)]) we have E A (a) - E A (rj) = lim /3 ^ 00 [i? A (cr) - 
H\(r])]/f3 = lmig^oofT//?) log[PA(cr, r])/P\(rj, a)]. Now, the fact that a and 77 form a stable pair 
implies limp^^ P A (a,rj) = lim^oo P\(r], a) = 1; hence E\(a) = E\(rj). By using results in Table 
[T]one can show that H\(a) and H\{rf) differ for a quantity exponentially small in f3. 

The remark above and the detailed balance condition suggests that the system cannot be 
trapped in loops longer than two. Indeed, consider a sequence a\,...,a n 6 S such that a,i + \ = 
To",, for all i = l,...,n — 1, and suppose, by absurdity, that o\ = Ta n . The property above 
implies that either E\(a{) — E^{a n ) = c > or E\{a\) — E\(a n ) = 0. In the first case from the 
detailed balance we get \ \og[P\{o-i,a n ) / P\{a n ,a\)\ — cf3\ — > in the limit (3 — > 00; hence using 
the hypothesis o~\ = To~ n we easily get an absurd. In the second case the detailed balance implies 
I log [Pa , o- n )/Px(a n , <7i)]| — > 0, that, togheter with o\ = Ta n , gives P\{ai^a n ) — > 1, which is 
absurd because by hypothesis we have 02 = T01. 

We say that a G S is a trap if either a is a stable configuration or the pair (a, Ta) is a stable 
pair. We also let M. C S the collection of all the traps. Now, we give a full description of the 
stable pairs in <S+i U Sc U S-\ (see Fig. 13. 2\ : the most general stable pair living in a sea of minus 
is made of rectangular flip-flopping droplets of chessboard plunged in the sea of minuses and well 
separated stable droplets of pluses living inside the sea of minuses or inside a chessboard droplet. 

Proposition 3.3 i) For any a G <S+i \ {+1} the pair (a, Ta) is not a stable pair, ii) Given C G C 
and a G Sc the pair (a, Ta) is a stable pair iff there exist k > pairwise non-interacting rectangles 
Re 1>mi , Re 2 ,m 2 , R i k ,m k , such that 2 < £ { < mi < L - 2 for any i = l,...,k, cjsr = +1^ 
(a coincides with +1 inside the rectangles) and o"a\sr = Ca\sr (°~ coincides with the chessboard C 
outside the rectangles), where 5? := (jf=i R^m- Given a G S-\ the pair (a,Ta) is a stable 
pair iff there exist k > 1 rectangles -Rf limi , i?^ 2jm2 , . .., R^ kiTnk! with 2 < ^ < nii < L — 2 for 
any i = l,...,k, and there exists an integer s G {1, . . . , k} such that the following conditions are 
fulfilled: 

1. Rt^rm n Re^mj = and U>2 for any i, j G {1, . . . , k}; 

2. for any j G {l,...,s} the family {Re^mj , Re a+1 ,m 3+1 , Rt k ,m k } is a family of pairwise non- 
interacting rectangles; 

3. o"a\sr = — Ia\sr where 3ft := Uf=i Rti,mi (°~ coincides with —1 outside the rectangles); 
4- a R f = +Ir« f or any j £ {s + 1, . . . , k} (a is plus inside Re a+ll m a+1 ,. ■ ■ ,R^ m J; 

3 ' rn 3 j ' rn 3 

5. for any j G {1, . . . , s} there exist k' = k'(j) > rectangles R' e ^ m ^ = ^/ )m / (j), R e' ,, m > , = 
Rip m , (j) such that the following conditions are fulfilled: 

5-1- R'^mt c R ij,mj for any i G {1, ... , k'}; 

5.2. for any j = 1, . . . , s the family {Rl , : i=l,...,k'} (recall R', , = R', , (j) for any 

i ' i i ' i i ' i 

i = 1, . . . , k' = k'(j) ) is a family of pairwise non-interacting rectangles; 
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A 





A 



Fig. 3.2: On the left (resp. right) the most general a E S^i (resp. Sc) such that (cr, Tit) is 

a stable pair. 



5.3. aw = where W = := (jf=i ^.mj 

54. either a^ m ^ = ("; ; ^ R , or ^ ^ = C^ ^,; 

6. for any i,j £ {1, ...,s} the two rectangles Rt-, m j an d ^.m* 'must be non-interacting if 
3.3. Basic tools 

In this Section we discuss the main tools that will be used in the following: first of all we notice that 
in our model the difference of energy between two configurations a,rj € S is not sufficient to say 
if the system prefers to jump from a to 77 or vice versa. Indeed, there exist pairs of configurations 
cr, 77 £ S such that the system sees a sort of energetic barrier both in the a — > rj and in the rj — > a 
transition. Let us define a sort of "communicating height" H(a, rj) for each pair (a, 77) £ S x S of 
the configuration space such that 

P A (<T,rj) =: e-I^i)-^)] . (3.3) 



More precisely, we consider a new hamiltonian H : S x 5 U S — > M defined as in (|2.4[) for any 
a G S and such that 

H(a,n) :=H(a)- log P A (a,rj) . (3.4) 

Note that, by virtue of the detailed balance principle (12. 5|) . we have H(a,rj) = H(rj,a). Remark: 
if either Pa(<t, rj) or -Pa(t/> cr) tends to zero in the limit (3 — ► 00 then i?(cr, 77) = max{if(<r), H(r))} + 
o(exp{— /3c}), for some strictly positive constant c; in other words in these cases the energetic 
barrier seen by the system is exactly the difference of energy between the two configurations. 

We notice that in [OS] it has already been remarked that the communicating heights allow to 
define the most general kind of reversible dynamics (see Section 3 in [OS] ) . Now we want to restate 
in this setup some of the results of [QSj that will be our basic tools in next sections. 

We say that a configuration a £ S is a local minimum of the energy iff H(a, rj) — H{a) > 
for any rj £ S \ {cr}. The local minima of the energy are nothing but the stable configurations 
defined above. A sequence of configurations uj = {u>o, . . . ,w n } is called a "path"; is the number 
of configurations in the path. We call "height along the path u" the real number 

& w := max H(ui-i,cOi) . (3.5) 
i=l,...,|w| 
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Given two configurations a,r] G S we denote by 0(tr, 77) the set of all the paths u = {ujq, . . . ,uj n } 
such that ujq = a and u n = 77. The "minimal height" (minmax) between a and 77 is defined as 

$((7,77) := min = min max H(uji-i,uJi) . (3-6) 

^G8(o-,rj) oje©(cr,?7) i=l,...,|o;| 

We remark that the function <E> : S x S — ► K. is symmetric, namely <I>(cr, 77) = $(77,0") for any 
a,rj € S. 

We give, now, the important notion of cycle: we say that A C S is a cycle iff for each cr, 77 G .4 

$0, 77) < min 4>(cr,C) • (3.7) 

In other words starting from any configuration in the cycle .4, the energetic barrier that must be 
bypassed to visit any other configuration in A is smaller than the one seen to exit the cycle itself. 
Given a cycle A C S we denote by F(A) the set of the minima of the energy in A, namely 

F(A) := {a G A : min #(77) = H(a)} ; (3.8) 

rjeA 

we also write H(F(A)) = -£^(77) with 77 G F(A). Given 77 G F(A) we define 

$(A) := min $(??,C) (3.9) 
(eS\A 

(it is trivial that 3>(«4) does not depend on the choice of 77 G F(A)), and the set 

U{A) :={(£S\A: 3r]€A such that H{rj, () = <S>(A)} . (3.10) 

Now, for any 7/ G 5 let be the probability over the process when the system is prepared in 
do = 77 and 

t v : = inf{n > : a n G T>} (3.11) 
for any V C S. We restate, without proof, some of the results of [OS] that we will use in the sequel: 



Lemma 3.4 Given Q C S, let a G Q and a' G S \ Q such that: i) there exist a* G S \ Q and a 
path u = {ujq = o~,...,u n = a*} such that u>i G Q and H(uJi^i,uJi) < H(uj n ^i,uJ n ) =: T for any 
i = 1, . . . , n — 1; ii) there exists a path ui' = {u' = a*,... ,u' n = a'} such that u)\ G S \ Q and 
H(u>' i _ 1 ,tL>' i ) < T for any i = 1, . . . , n; Hi) ming-g^gs^ H (er, 77) > V if and only if a = u) n -\ and 
n = u> n . If we define 

A := {77 G S : 3u> = {loo = 77, . . . , u n = a} such that U)\,..., ui n -i G Q and <& w < T} 

then i) A C Q; ii) A is a cycle with &(A) = T and a* G U(A); Hi) &(o~, a') = T (that is T is the 
minmax between a and a'). 

Lemma 3.5 Given a cycle A C S, 

i) for all e > and for all a G -^(-4) 

P CT (exp{d>(^) - H(F(A)) - (3e) < t s \ A < exp{<S>(A) - H(F(A)) + /3e}) ^ 1 (3.12) 

ii) there exists 5 > such that for any a,T] G A 

IV (t„ < TS \ A , r n < exp{$(„4) - H(F(A)) - (35}) ^ 1 (3.13) 
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Hi) for any a £ A 

P CT (v Ts ^ G U(A)) ^ 1 (3.14) 
iv) for any a £ A, rj £ U(A), e > and /3 sufficiently large 

^ (°r SXA = v)> (3.15) 

We note that the lower bound on ts\_a_ in the statement i) in Lemma [3.51 is an easy consequence of 
the reversibility property: 

Lemma 3.6 For any a, n £ S such that <3?(<7, rf) — H(a) > and for any 5 > 

P CT (r„ > exp{$(cr, n) - H(a) - (35}) ^ 1 . (3.16) 

Few important remarks which are very peculiar of our PCA model. As it has been noticed 
above if the system is in the state rj £ S, then there exists a unique configuration where it jumps 
with high probability. This configuration has been denoted by Try. Thus, given rj G S we define the 
downhill path starting from n as the unique path lo = {loq, . . . ,ui n } such that ojq = 77, Tu^i = uji 
for any i = 1, ... ,n, and oj n is a trap; we also set rj := u n . We remark that to each n G 5 we 
can associate either a unique stable configuration or a unique stable pair. We define the basin of 
attraction of a trap r/ G A4 as the set 

Bfo) := {C G 5 : C = "} (3.17) 

and the truncated basin of attraction £>(<r) C 13(a) as the set of all the configurations r\ G 13(a) such 
that 

$(ri,a) < min $(77, • (3.18) 

(£S\B( V ) 

It can be easily proven that B(rj) is a cycle. 

In the following we will often have to evaluate T(r?) := min^g^yg^) <3?(n, Q for some trap 77 G M.. 
A convenient way to proceed is the following: say that a path uj = {ujo, . . . , uj n } is uphill iff the path 
uj 1 = {cjq, . . . , ui' n }, where oj[ = oj n -i for any i = 0, . . . , n, is downhill. Consider the set £(77) of paths 
uj = {ojq, . . . , w n } such that u = 77, {ojq, . . . , w n _i} is an uphill path in B(rj) and u> n £ S \ B(n). 
We remark that T(n), in words the barrier that must be bypassed to exit from the basin B(rj), is 
given by 

1(77) = min $ w . (3.19) 

3.4- Behavior of traps 

In this subsection we clarify the geometrical conditions for the shrinking or the growing of a trap, 
that is we study the evolution of the system prepared in a stable configuration or in a stable pair. 
We let A := [2/h] + 1. 

We first consider the case of a single rectangular droplet of chessboard or pluses in the sea of 
minuses; we show that if the droplet is small enough, namely its shortest side is smaller than A, 
then it tends to shrink, otherwise it tends to grow. 

Proposition 3.7 Let ( £ {+1, C e ,C°}; 77 G M. such that there exists a rectangle Ri )T n, with 
2 < t < m, such that n A \-n = — l A \-p and n-n = Cs ■ Thus 
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Fig. 3.3: Growth of a chessboard droplet inside the sea of minuses: appearing of a 

protuberance. 



i) if £ < X, then n is subcritical, that is P^(r_i < t$ C [js +1 ) — ► 1, and for any e > 

P^ (exp {2(3h{£ - 1) - (3e} < r_i < exp {2f3h{£ - 1) + 0s}) ^ 1 ; (3.20) 

ii) if £ > X, then r] is supercritical, that is W v {i~s c uS + i < - — 1, o,nd for any e > 

P„ (exp {2/3(2 -h)- fie] < t Sc us +1 < exp {2/3(2 - h) + fie}) ^ 1 . (3.21) 

Similar results can be stated in the case of a single droplet trap plunged inside the sea of 
chessboard. 

Proposition 3.8 Let C G C and r/ G M a trap such that there exists a rectangle Ri m , with 
2 < £ < m, and r7 AV -„ = C AV -n and rps = +l-p . Thus 

i) if £ < X, then P v (tc < r+i) \ ) an d f or an y £ > 

P,, (exp {2f3h(£ - 1) - /3e} < r c < exp {2j3h{£ - 1) + /3e}) 1 ; (3.22) 

ii) if £ > X, then P^(r+i < tq) - — ? 1, and /or any e > 

(exp {2/3(2 -h)- I3e) < r +i < exp {2/3(2 - /i) + fie}) 1 . (3.23) 

Now, we give two heuristic arguments supporting the Propositions above. We consider the case 
£ = m even and 77^ ( = : by using ()2.9p one can show that the energy of rj, with respect to 

the configuration —1, is 

E\(rf) - E h k (-l) = -2h£ 2 + 8£ . (3.24) 

Thus, E^(r}) — E^(— 1) is a parabola whose maximum is achieved &t £ = 2/h suggesting the 
conjecture that the critical length is 2/h. 

A dynamical argument strengthens this conjecture. Consider the most efficient growing mech- 
anism: from results in Table [1] this mechanism is the appearance of a single plus protuberance 
adjacent to one of the four sides of the rectangle. The probability associated to such an event is 
exp{— 2/3(2 — h)}, so that the typical time to see this event is r gr ~ exp{2/3(2 — h)}. In Fig. 13.31 it 
is shown that once the protuberance has appeared on one the four sides of the rectangle, with high 
probability a new slice is filled with chessboard. 

Now we have to look for the efficient shrinking mechanism: Table [T] suggests this mechanism is 
the "minus corner persistence" , that is a minus spin on one of the four corners of the rectangle is 
kept fixed during the flip-flop of the stable pair. By reiterating this mechanism £ — 1 times (see 
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Fig. 3.4: Shrinking of a chessboard droplet inside the sea of minuses: persistence of a 

minus corner. 



Fig. I3.4|) a full slice of the droplet is erased. In terms of probability each step costs exp{— 2(3h}, 
hence the typical shrinking time is r s h ~ exp{2(3h(£ — 1)}. By comparing r gr and r s h we find that 
the critical length should be 2/h. 

A similar argument can be done in the case of a plus droplet: the growing mechanism is still 
the formation of a protuberance. About the shrinking mechanism: after a first step of "corner 
erosion", like in the standard Glauber case a corner is flipped into minus, one minus spin appears 
at the corner. The best thing to do, as a second step, is to flip simultaneously both the minus spin 
and its adjacent plus spin. This event costs still exp{— 2/3h} and results in a shift of the minus 
"lacuna" on the side of the rectangle. By iterating this mechanisms a sort of merlon is formed and 
a stable pair is reached in a typical time r s h ~ exp{2(3h(£ — 1)}. 

A stronger version of the above Propositions can be proved; it is possible to describe in detail 
the way in which droplets shrink or grow. Indeed we state: 

Proposition 3.9 i) Let C E C, n E Ai such that there exists a rectangle Re, m , with £ < m, such 
that VaVr^ = —lA\R e an d n R e = ^R e ■ Let A 1 the set of traps a E M such that there exists 
a rectangle R with side lengths (£,m + 1) or (£ + 1, m) such that R D Ri t m, °~a\r = ~—A\R an< ^ 
a-ft = C-ft. Let A" the set of traps a E M such that there exists a rectangle R with side lengths 
(£, m — 1) such that R C Ri^ m , ° a\r = ~-A\R an( ^ a ~R = V ^ — ^ then 

T(n) =H(n) + 2(3(2 - h) and P„ ( a TsXf3{?)) E (J B(a)) ^ 1 , 

that is starting from r\ the system exits B(n) and enters into one of the basins 13(a), for some 
a £ A'. If £ < X then 

T(n)=H(n) + 2(3h(£-l) and P,( %5h) e/) ^ 1 , 

that is starting from n the system exits B(rj) and reaches directly one of the traps a in A" . 
ii) Let C E C, r/ E M such that there exists a rectangle i?^ m , with £ < m, such that VA\R e = 
^A\R t an d r fR t , = ~^—Ri ■ ^t ^ ^ e se ^ °f traps o~ £ M such that there exists a rectangle R 
with side lengths (£, m + 1) or (£+l,m) such that R D Re, m , &a\r = ^aVr an< ^ a ~R. = ^-R- Let A" 
the set of traps a E M such that there exists a rectangle R with side lengths (£, m — 1) such that 
R C Re, m , CA\fl = ^a\r, an d a R = +-R- V X < £ then 

T(r l ) = H(r / ) + 2(3(2 -h) and P, %%) G (J B(a)) ^ 1 , 

V aeA' ) 
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that is starting from rj the system exits B(rj) and enters into one of the basins B(a), for some 
a £ A'. If £ < A then 

T(rj) = H(n) + 2(3h(£-l) and P„ (a Ts ^ } G A") 1 , 

that is starting from n the system exits B(rj) and reaches directly one of the traps a in A" . 

Finally, we state under which conditions a general trap shrinks. A trap is made of rectangles of 
pluses or chessboard inside a minus or chessboard sea. The idea is that the configuration shrinks 
iff each single rectangular cluster shrinks. Note that the following proposition strictly contains 
Propositions 13.71 and 13.81 

Proposition 3.10 i) Let C G C and n G Mc- There exist k > 1 pairwise non-interacting rectan- 
gles Rt 1)mi , Re k ,m k such that 2 < li < m < L - 2 for any i = 1, . . . , k, rj$t = +1 R and n A \ n = 

Ca\5R where 51i := |Ji=i Rli,m % - Thus, if li < A for any i = 1, . . . , k, then F v (tc < T+l) — ^ 1> an d 
for any e > 

P,, (exp {2(3h(£ - 1) - 0e} < t c < exp {2(3h(£ - 1) + fie}) ^? 1 , (3.25) 

where £ := max{4, • • • , £k}- If there exists j G {1, . . . , k} such that £j > A, then P 7? (r + i < tq) - — 
1, and for any e > 

P^ (exp {2/3(2 - h) - {3e} < r +1 < exp {2/3(2 - h) + f3e}) ^ 1 (3.26) 

(note that in the case k = 1 we recover Proposition \3.8\) . ii) We consider, now, a situation where 
in a sea of minuses there are rectangles of chessboard and, possibly, rectangles of pluses inside the 
sea of minuses or inside the chessboard droplets. More precisely, let n G M. We suppose that there 
exist k > 1 rectangles Ri lt7ni , ■ ■ ■ , Rt k)TrLk , with 2 < £i < rrii < L — 2 for any i = 1, . . . , k, and 
there exists an integer s G {1, . . . , k} such that the conditions of point (Hi) in Proposition \3.3\ are 
satisfied (note that in the case s = the trap is a stable configuration; in the case s = and k = 1 

Proposition \3. 7| is recovered). Thus, if £-i < A for any i = 1, . . . , k, then P r/ (r_i < ts cU s +1 ) - — ^ 1, 
and for any e > 

P^ (exp {2f3h{£ - 1) - fie} < r_i < exp {2f3h(£ - 1) + /3e}) 1 , (3.27) 

where £ := max{4, • • • , £k}- If there exists j G {1, . . . , k} such that £j > A, then n is supercritical, 
that is H s v {ts\s_ 1 < T-i) - — ? 1, and for any e > 

P r , (exp {2/3(2 - h) - /3e} < r SV ;_, < exp {2/3(2 - h) + fSs}) 1 . (3.28) 

3.5. Exit from the metastable phase 

We can now give the theorem describing the exit from the metastable state. Suppose that the 
system is prepared in the metastable state, do = —1, in the following theorem we state that the 
first exit time r+i is exponentially large in /3 and we find its order of magnitude. Moreover, we 
state that before reaching +1 the system visits C = {C°, C e } and that the typical time to jump 
from — 1 to the chessboards is the same as the time needed to jump from the chessboards to +1. 
More precisely, let us denote by Q~i C S the set of configurations n G S such that there exists a 
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rectangle R\ t \ such that ??a\^ a x = — Ia\;r a a an< ^ 7 /# A A = a with C £ C. Let 7? G Q-i, we call 
"protocritical height" the energy 

r := ££(77) - £a(-1) + 2/i(A - 1) = -2hX 2 + 2A(4 + /i) - 2/i , (3.29) 

where the second equality follows from (|3.24|) . In some sense j3T is the communication height 
between the largest subcritical droplet and the smallest supercritical droplet; more precisely: let 
C_i C S the set of configurations 77 G 5_i such that there exists a rectangle -Ra-i,a such that 

VASRy-^ = -IaNHa-x.a and = °Ha-i,a With C G C - 

Let 77 G O-i : we call protocritical droplet corresponding to 77 one of the configurations obtained 

by flipping in 77 a minus spin external to R\—i x an d adjacent to one of the plus spins of the internal 

chessboard and all the spins associated to sites inside -Ra— i,a- We let 7r_i(77) the set of protocritical 

droplets corresponding to 77 and V-\ '■= U^ g o_ 1 7r_i(7/), the collection of the protocritical droplets. 

It is easy to check that H(rj, £) = T for any 77 G 0-\ and C, G vr_i (77). 

Theorem 3.11 With the notation introduced above: i) the system visits V—\ before visiting C, 
namely 

P-ifa-i < r c ) ^ 1 ; 
ii) the system visits C before visiting +1, namely 

P_i(r c < r+i) ^ 1 ; 

Hi) for any e > 

P_ l(e /3r-/3 £ <TC< e /3r +/3£) /s^oo x . 

7tj) /or any e > 

P_i(e^ £ < r+i < eF^) ^ 1 . 

The proof of Theorem 13.111 will be the argument of Section HJ 

In the above theorem we have stated that during the exit from the metastable — 1 state, the 
system visits the competing metastable state C and, finally, reaches the stable state +1. Now, we 
want to give a more precise description of the path followed by the system during its exit from — 1; 
first of all we define a suitable tube T-\. Let 0_l := Q-\ the set of A x (A — 1) chessboard droplets 

in the sea of minuses. For each k = 0, 1, . . . , 2A — 6 we define recursively the sets 0_\ (j/* - ^), 
where rj^ k ~ l ^> G 0^ 1 \r/^ k ~ 2 ^): let rf k ~ x ^ G 0^i(rf^ k ~ 2 ^) a configuration such that there exists a 
rectangle Rp m , with 2 < £ < m, such that 77^ ^ = C-s , with C G C, and 77.^ — = — l A \"p 
Then we define 0_l(j]^ k ^) as the collection of configurations £ such that there exists a rectangle 
Re, m -i such that Rt, m -\ C Re, m , Ca\r, , = ~Ia\r, 1 and Cr, , = c r, , > with We 

j i\. \iif. m - 1 yi^m— 1 ll £ f m — 1 ll £,m— 1 

remark that 0^ -7) (^( 2A ~ 8 )) is a set of 2 x 2 chessboard droplets and of 2 ^ 6 ^(r/ 2A 7 )) is made of 
two configurations with a plus spin in the sea of minuses. 

We note that Proposition 13.91 implies that the process enters, with high probability in the 
limit j3 — > 00, into the set 0_i (77®), when it exits from the basin of attraction B(rji) with 

m g oSfaf*- 1 )). 

Now, given the 2A — 5 configurations rji G for any £ = 0,1,..., 2A — 6 and the 

2 A — 6 integer numbers *!<•••< t2A-6 we se t *o — an d *2A-5 — ^2A-6 + 1> an d we say that 
a path u = {u to ,...,uJt 2X _ B } belongs to the set T (77W, . . . , 7?( 2A-6 ); t u . . . , t 2 \-s) iff oo tl = V® 
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and oj tl , ■ ■ ■ ,u)t i+1 -\ & B(r/^) for any i = 0, . . . ,2A — 6 and W( 2A 5 = — 1. Note that a path in 
7" (rf°\ . . . , ry( 2A ~ 6 ); ti, . . . , t2x-6) is one of the "standard" shrinking paths that the system follows 
when a chessboard droplet r/°) G shrinks. More precisely, given r/ ) G cf°i, we define the 
tube 



2A-6 



V : =U U U ^(^ (0) ,---^ (2A - 6) *2A-e) • (3.30) 

i=l jjWgOW^Ci- 1 )) tl<— <t 2 \-6 

In other words 7 l(o) is defined as the set of paths =G S(t/°), —1) such that there exist 2A — 6 
configurations 77 W G 0_{(r7^ — ^) with i = 1, . . . ,2A — 6 and 2A — 6 integer numbers <]_<•••< i2A-6 
such that uj G 7" (r/°) , . . . , r/ 2A_6 ) ; ti , . . . , *2A-6) • We state, now, the following Lemma: 

Lemma 3.12 Let r/ G 0-\, we have 



■ n(0) 



i/te trajectory {ctq, . . . , er T _i } G 71 ( ) ) ^ — ^ 1 



Proof of Lemma \3.1Si The Lemma easily follows by applying recursively the Proposition 13.91 and 
the Markov property. □ 

Finally we define the exit tube TLi: a path to = {ujq, . . . ,uj n } is an element of 71_i iff there 
exist 2A — 5 configurations 77W G ©^(t/ 4-1 )) with i = 0, 1, . . . , 2A — 6, the integer numbers t\ < 
■ ■ ■ < i2A-6 = n — 1 and a path u/ = {uj' , . . . ,Lo' n } G T (rj(°\ . . . , ?^( 2A_6 ); t , . . . , *2A-6) such that 
uji = u)' n _i for any i = 0, . . . , n. In other words 71. 1 is the set of paths obtained by time reversing 
one of the standard shrinking paths associated to the droplets in 0-\. 

Theorem 3.13 Let at be the process started at —1, let f_i := max{£ < t s \^_ 1 : at = —I}, then 

P-i ( cr r S \ A _ 1 e the trajectory |o- f _ 1; cr f _ i+ i, . . . , cj 7 - sv ^_ i _i| G T-A 1 . 

Proof of Theorem \3.13l The theorem is a straightforward consequence of the time-reversing argu- 
ment (see \OS\ E] ) and Lemma 13.121 □ 



4. The minmax between the metastable and the stable state 

The proof of the theorems describing the exit of the system from the metastable state is based 
on the general lemmata given in the Subsection \3. 31 The highly not trivial model dependent part 
consists in finding the minmax between the metastable and the stable phases. It is clear that new 
ideas must be used to answer this question in the case of a parallel dynamics with respect to the 
Glauber case. Indeed, the fact that the system can jump from any configuration to any other, 
highly complicates the structure of the possible trajectories in the configuration space. 

First of all we define a sort of generalized basin of attraction of — 1: let C <S_i the set 

Q-i := {a G S-i : a = — 1 or a subcritical} , (4-1) 

where a subcritical means that a is a trap such that there exist k > 1 rectangles on the dual lattice 
i?< limi , . . . , i?£ fejmfe , with 2 < $4 < nii < L — 2 for any i = 1, . . . , k, there exists s€ {0,...,fc} such 
that the conditions of point in) in Proposition 13.31 are satisfied (note that in the case s = the 
trap is a stable configuration) and ti < A for any i = 1, . . . , k. To fix the ideas: if a consisted of a 
single chessboard rectangle in the sea of minuses, then its shortest side length t should be smaller 
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than A. The set Q_\ is a sort of "generalized" basin of attraction of the state — 1_, in the sense that 
for any a G Q-\ the process started at a would visit — 1 before exiting with high probability 
in the zero temperature limit, namely 

^ (r-i < r s \g_^ ^ 1 • (4.2) 

In the following Lemma we state the main properties of the basin Q-\\ Y is the minimal energy 
barrier that must be bypassed to exit Q-\\ a minimal exit path from Q_\ reaches S \ Q-\ in a 
protocritical droplet rj G V-\- 

Lemma 4.1 Let r/ G V—\, i) there exists a path uj = {luq = —1, ... ,u> n = r]} such that uji G Q~\ 
and H(uJi^i,uii) < H(uj n ^i,uj n ) = H( — l) + j3T for any i = 1, . . . , n — 1; it) there exists a path 
uj' = {uj'q = r],...,Lo' n G C} such that G S \ Q~\ and H(uJi-\,u)i) < H(—l) + f3T for any 
i = l,...,n. in) $(C7_i) = H(-l)+(3T; iv) for alia G Q-\ and rj G S\C?-i, H(a,rj) = 0T+H(-1) 
if and only if a G C_i and r} G -K-\[a). 

We postpone the proof of the above lemma to the end of this section. Let us define the set 

A-\ := {rj G S : 3to = {loo = rj, . . . ,ui n = —1} such that 

cj u . . . , w n _i G Q-i and $ w < H(-l) + /3T} 

From Lemma [3.41 and Lemma |4. II we have that A-\ is a cycle, A~i C G-i, &(A-\) = H(—l) + f3T, 
U(A-i) D V-x and $(-l,C) = /3r + #(-1), with C G C (that is /?r"+ il(-l) is the minmax 
between — 1 and C). 

Proof of Theorem \3.11[ Let be the process started at — 1. We first try to describe the exit from 
the basin by applying Lemma 13.51 and using recurrence in G-i- We firstly remark that from 
the reversibility Lemma we have that for each e > 



P-i [rs\ S ^ > e") ^ 1 (4.3) 

Now, from item iv) in Lemma 13.51 from the definition of Q-\ and from Proposition 13.91 we have 
that for each a G Q-\ and each e, 5 > 

P CT (at < e^ +f3 \ a t eS\ G„x) > (4.4) 



-OC 



For any e > 0, we set T(e) := exp{/?r + /3e}, N{e) = [exp{/3e/2}] — 1 G N, and consider the intervals 
I k (e) ■= T(e)exp{-/3e/2}[/c,/c + 1) for any k = 0, . . . ,N(s). Then 



-I [rsXS-, < e W ) = l-P-i(r 5W _ 1>e ^ 

JV( £ ) 



fe=o (4.5) 
1_ sup [l -P,(3t< e ^ £ / 2 , a t GcS\g_ ^ 



> 1 



N(e) 



1 - e-^/ 2 

where we have used ()4.4p . From (|4.5p we get the upper bound on the exit time 

P-i (r S \g^ < e^) ^ 1 (4.6) 
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for any e > 0. Now, by using the reversibility Lemma we have that 



-l 



e P-i) ^ i (4.7) 



and recalling (|4.3p and (|4.6p we ca state that the system prepared in — 1 exits <5_i through V-\ 
in a typical time exp{/3r}. Moreover, C C S \ Q-\ (the chessboards do not belong to G-i) and 
equation (|4.7p imply that P-i is visited before C, more precisely 

P-i {tr_ x < re) ^ 1 , 

completing the proof of the statement i) in Theorem 13.111 

Now, we use the Markov property to restart the system in some configuration of V-i- Point i) 
in Proposition 13.91 directly implies point ii) in Theorem 13. Ill Moreover, point Hi) in Theorem 13. Ill 
is easily proven by remarking that [3T > 2/3(2 — h) and by using Proposition 13.91 

Up to now we have described the jump from — 1 to the chessboards. Now we use the Markov 
property to restart the system in C G C and we prove point iv) in Theorem 13.111 by following the 
same scheme used above. We just sketch the proof: let Qc C Sc be the set 

Qc '■= {c G Sq '■ d G C or a subcritical} , (4.8) 

where a subcritical means that a is a trap such that k > 1 pairwise non-interacting rectangles 
Re 1>mi , • • • , Ri k ,m k such that 2 < ii < mi < L — 2 and £i < X for any i = 1, . . . , k, ijgi = +1sr and 
^AXK = Ca\R where 3ft := Ui=i ^tum an d C G C. To fix the ideas: if a consisted of a single plus 
rectangle in the sea of chessboard, then its shortest side length I should be smaller than A. Then 
we state the analogous of Lemma 14.11 let us denote by Qc the set of configurations rj G S such 
that there exists a rectangle R\ t \ such that Va\r x a = ^a\r x x and ^Rx x = ~^~-R\ x w ^ n C € C. 
We denote by Vc the set of configurations r\ G S such that there exist a rectangle R\-i t \, a site 
x G A \ ^, adjacent to one of the two sides of R\-i,x of length A, and C G C such that 
C(x) = -1, Va\(r XiX u{x}) = C A\(R x , x u{x})i ^ A , A = +Iii A , A and r ?( x ) = + L 

Lemma 4.2 Xei C G C, r/ G Pe? i) ^ere exists a path uj = {ujq = C,Ui, . .. ,co n = rf\ such that 
&i G Qc and H(uji-i,uJi) < H(u} n -i } u) n ) = H(C) + f3T for any i = 1, . . . , n — 1; ii) there exists a 
path to' = {cu'q = 7/,. . . ,(jj' n = +1} such that u[ G S\Qc and H(ui-i,u>i) < H(C) + (3T for any 
i = l,...,n. Hi) §{G C )=H{C) + PT. 

As before the statement iv) in Theorem 13.111 follows from the Lemmata 14.21 13.41 13.51 and point ii) 
in Proposition 13.91 □ 



Proof of Lemma \4-1\ We start by proving point i): let us consider a protocritical droplet rj G V—\ 
and the Ax (A — 1) chessboard droplet 7/ ) G 0—\ such that 7t_i(t/ )) = 77. First of all we note 
that: H(r),r)W) = H{-l)+(3T. 

Now, recall O^l = 0_i and consider a sequence of configurations r/ 1 ) G 0_l(rj^),. . . , r/2A-6 G 

of 2 1 A ~ 6) (r ? ( 2A - 7 )). From Proposition [3Jj] we have that for any i = 1, . . . , 2A — 7 the barrier to exit 
the related basin of attraction is 

T(t/W) < H(ti®) + 2l3h{\ - 2) < H{rf ^) + 2ph(X - 2) 

< H(i]^) + 2/3(2 -h) = H(-l) + /?r 
Note that for i = 0, 1 the first inequality is, indeed, an equality. 
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The above inequalities allow to construct a path {ojq = rj(°> , . . . ,uj n = — 1} G 7^( ) connecting 
r/°) to —1 and such that H(oJi,uji + 1) < H(— 1) + /3r for any i = 0, . . . , n — 1. Finally, we remark 
that the path {wq = u; n ,u4 = w„_i, . . . ,ui' n = ujQ,oj' n+1 = rj} satisfies the properties of point i). A 
similar construction can be repeated for point ii). 

Now, we come to the main points Hi) and iv): our goal is to prove that &(G-i) = H(— 1) + f3T. 
First of all we notice that point i) above implies 

$(G-i) < pr + H(-i) , 

hence our calculation is reduced to prove the lower bound &(G-i) > f3T + H{— 1). We have to 
examine all the paths connecting G-i with S \ Q—\: such a path {loq,uji, . . . ,uj n } has at least a 
direct jump from G~i to S \ G-i, that is there exists k G {0, . . . , n — 1} such that u>k G G-i and 
uJk+i G S \ G-i- Hence, for any ui connecting G-i with its exterior we have 

Qui > nain H (a, r]) . 

&eG-i,ri£S\g-i 

Given a G G-i, we consider the configuration rj G S\G-i that can be reached with the smallest 
energetic cost, namely rj G S \ G~i is such that 

H(a, n) = min H(a, C) (4.9) 

We have to prove that H(a,rf) > (5T + H(— 1) with the equality valid if and only if a G C_i and 
7] G 7r_i(cr). 

First of all we note that P\(a, if) - — ? otherwise we would have r\ G B(a). From Table 1 we 
get that there exists x G A such that 

logp x (7/(x)|(7) = -2/3(2 - /i) + o(e-^ c ) (4.10) 



for some positive constant c. We consider, then, £ = rf , and we remark that (14. 9j) implies £ G G-i- 
Note that 

log P A (^r?) = log P A (a, C) - logp x (C(x)k) +log Pa: (7/(x)|a) (4.11) 
that, together with (|3.4p . implies 

H(a,rj) = H (a, 0- log P A (a, V )+ log P A (a,Q 

= H(a, C) - log p x (r](x)\a) + logp a; (C(x)|cr) 

> H(C) - logp x ( V (x)\a) + logp x (((x)\a) (4.12) 

We can characterize £ as follows: by using Propositions 13.31 and 13.101 we have that there exist 
k > 1 rectangles P^ ljmi , . . . ,Ri kt m k satisfying the conditions of point (Hi) in Proposition 13.31 with 
respect to the configuration £, and such that li < \ for any i = 1, . . . , k. 

Let us consider C' G S such that C' k -j- = — 1, -5 and C' A \, lfc -75 = Ca\i ifc r 

S S Uf =1 %., m . - U i=l- R « i ,m i S A\U^ =1 i?^, m . S A\U i=1 ^, mi 

By recalling that there exist an uphill path joining £ to £ it is easy to prove that there exist s > 
rectangles R l m , . . . , Pi, m , with m k+i > 1 for all i = 1, . . . , s, such that £' „ = 

"\u i=1 «i im( , +i 

— 1 A \, |S and C' coincides with a chessboard or +1 inside Pi m , , . for all i = 1, . . . , s. 

— il \ u i = l-' t l,m fe + i — ' R + l 

We let 6° (resp. bV) the horizontal (resp. vertical) side length of the rectangle P^ iimi for any 
i = 1, . . . , k + s. We set b° := X)i=i &i > ^ := Si=i anc ^ we remark that 

+ = 6° + b v . (4.13) 
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We suppose, now, b° + b v > 2A. By a direct evaluation of the energy of the rectangles it is easy 
to show the bound 



E(()-E(-l) 



E(() - E(-l) 



+ 



> 



k+s 



E(C) ~ E(-l) + [E(C') - E(-l)] > 4(6° + 6") -2hJ2 W (4-14) 



i=l 



Where we have used that two of the -Ri, mi , with i = k + 1, . . . , k + s, rectangles can interact iff they 
are filled with different parity chessboards. Now, by using the subcriticality of £, namely by using 
li < A for any i = 1, . . . , k (recall l% = \ for all i = k + 1, . . . , k + s), we have 



k+s 



j2bZbv<(\-i)[b°+b v -(\-i) 



(4.15) 



indeed, let £ = maxj =lv . ^ £i, we have 

k+s k+s k+s 



b°b\ = <iJ2 m - + b v — £] < (X — l)[b° + b v - (\ - 1) 



i=l 



i=l 



i=l 



where, in the last inequality, we have used b° + b v > 2 A. Now, recall A = [2/h] + l = 2/h + e for 
some e € (0, 1), then the inequality 



2 - h(\ - 1) 



h{l-£)>0 , 



(4.16) 



(jUD, (|4TT5|) and the hypothesis b° + b v > 2 A imply 



E(()-E(-l) > 4(6° + b v )-2h{\-l)[b° + b v -(A-l)] 

= (6° + b v ) [4 - 2h{\ - 1)] + 2h{\ - if 

> 2A[4 - 2h(X - 1)] + 2h(X - l) 2 

= 8A - 2h\ 2 + 2h 



(4.17) 



Now, by using (|4.1Uj) . (|4.12p and (|4.17p we get that for any 5 > there exist (3 large enough 
such that 

H(a, i]) - iT(-l) > 8/3A - 2/?/iA 2 + 4/3 - 5 (4.18) 

Finally, from the equation above it follows, by choosing 5 small enough, that H (a, rf) > f3T+H(—l). 

We come, now, to the case b° + b v < 2 A — 1. First of all we notice that £ and r\ differ for the value 
of a single spin, this implies that in r\ there is a single supercritical rectangle. More precisely, there 
exist k > 1 rectangles R e i m > , . . . , m i satisfying the conditions of point (Hi) in Proposition 13.31 

1 ? 1 fc/ ) fc/ 

with respect to the configuration rj, and such that £[ > A and l\ < A for any % = 2, . . . , k'. 

Let 1° (resp. f) the length of the horizontal (resp. vertical) side of the supercritical rectangle 
^j.m^- We note that b° + b v < 2 A — 1 implies that 1° + £ v is surely less than 4A; one can show, 
indeed, that 1° + £ v does not exceed 2A + 4. Under this condition it is easy to show the bound 



E(rf) - E(-l) > 4(1° + £ v ) - 2h£° 



(4.19) 



indeed the energy of rj can be bounded from below with the energy of an £° x £ v chessboard droplet 
in the sea of minuses. 
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Now, let C G {C e ,C°} such that rj(x) = C(x) and consider the collection of rectangles TZ C 
{Ri umi , ■ ■ ■ ,Re k , mk ,Ri,m k+1 , ■ ■ ■ ,Ri,m k+s }, such that for each R £ TZ either = +1r or the 
chessboard part of R coincides with C. By remarking that TZ is a collection of pairwise not 
interacting rectangles, we can find a positive integer A < 9 such that if we let 



V := R e[ , m { \ (J R and ^ := 



fy.mj \ U R 



Ren 



\V\ (4.20) 



ReTZ 

we get the lower bound 

E(rj) - E(-l) > [E(rf) - E{-1)] + 2h(N - A) > 4(f> + 1) - 2h£°t + 2h(N - A) (4.21) 

where, in the last inequality, we have used (I4.19p . We remark that N is a lower bound of the 
number of sites in Rp m i not belonging to any cluster of £ that will persist in rj. 

We consider, now, the geometrical projection of the rectangles R£ im ., with i = 1, . . . , k+s, onto 
one of the horizontal (resp. vertical) sides of Rg> m ' . Such a projection is a collection of, maybe 
not disjoint, segments; we denote with p° (resp. p v ) the length of the union of these segments. By 
definition we have p° < b° (resp. p v < b v ). We mention the following interesting bound on N: 

N > l°t - p°p v (4.22) 

Indeed, P sites of V are associated to each unit segment of the horizontal side of m i not 
belonging to the the projection of the rectangles Ri itnii . Moreover, p° (not already counted) sites 
of V are associated to each unit segment of the vertical side of Ry not belonging to the the 
projection of the rectangles Rg umi . Hence, 

N = \V\>(£°- p°)f + {f - p")p° = i°t - p°p v . 

Without loss of generality we can, now, suppose b° < b v . This implies b° < A — 1, indeed if it 
were, by absurdity, b° > A, then it would be b° + b v > 2A > 2A — 1. We distinguish among four 
different situations. 

Case 1. p° < A - 2. By inserting (14221) in (14211) we get 

E(rj) -E(-l) > 4{£° + F)- 2M°t + 2h{£°t - p°p v ) - 2hA 

> 4{£° + £ v ) - 2h(X - 2)£ v - 2/iA (4.23) 

where, in the last inequality, we use p° < A — 2 and p v < £ v . Now, recalling £ v > A, for h small 
enough we have £ v > A, hence 

E{rj) - E(-l) > A{£° + 1) - 2h(X - l)t = U° + 2t[2 - h(X - 1)] 

> 4A + 2A[2 - h(X - 1)] = 8A - 2/iA 2 + 2hX > T (4.24) 

where we have used (|4. 16p . 

Case 2. p° = X - 1 (recall p° < X - 1) and (£°,£ v ) / (A, A). As in Case 1 we get 

E(r}) - E(-l) > 4(£° + t) - 2hp°p" - 2hA > 4(f> + t) - 2h(X - 1)A - 2/iA (4.25) 
where in the last inequality we have used p° = A — 1 and p v < A, indeed 

p v < b v < 2A - 1 - b° < 2A - 1 - p° = X 
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Fig. 4.5: The three possible situations that must be taken into account in Case 4. 



Hence, by using 1° + £" > 2 A + 1 we get 

E(rj) - E(-l) > 4(2A + 1) - 2h{\ - 1)A - 2/iA 

= 8A - 2/iA 2 + 2h{\ - 1) + 2/i(A - A) + [4 - 2/i(A - 1)] 

= T + 2h(X - A) + [4 - 2h{\ - 1)] > T (4.26) 

where we have used 4 — 2/i(A — 1) > (see inequality (|4.16p and have chosen h small enough in 
order to get A > A. 

Case 3. p° = A — 1, t° = £ v = A and k + s > 2. In this case we have the easy estimate (sec [KO] 
and [NO]) N > 2A, hence from (|42TD we get 

E(rf) - E(-l) > 8A - 2h\ 2 + 2h{2\ - A) . (4.27) 

By choosing h small enough we get E{rf) > V + E(—l). 

Case 4- P° = A — 1, 1° = £ v = A and A; + s = 1. We have to consider the three configurations rji, r\i 
and 7/3 depicted in Fig. 14.51 where the A x (A — I) rectangle is filled with one of the two chessboards 
and r](x) = +1. By a direct evaluation of the energy we have that 

E(rj 2 ) = T-4h, E{rn) - E{r] 2 ) = 4 + 2/i > Ah and Efa) -£(%) = 2h(\ - 1) > Ah . 

Now, starting from 772 the lowest energy jump toward Q-\ consists in reverting all the spins inside 
the rectangle R and the plus spin associated to x. The cost of such a jump is exp (—A(3h), indeed we 
have to pay in order to keep the two minuses around x. Finally we have that H(a, 77) = f3T + H(—V) 
iff a G O-i and 77 G 7r_i(<r). □ 



5. Proof of the Propositions 

In this section we prove the Propositions stated throughout the paper: the tools which will be used 
are those outlined in Subsection 1 3. 3\ 

The Proposition 13.11 is a straightforward consequence of Lemma 13.21 and the characterization of the 
local minima of the energy given in [NSlj . 
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Proof of Proposition \3.3l We just give a sketch of the proof, i) Let a G <S + i \ {+!}■ There exists 
x G A with two neighboring sites in the sea of pluses such that o~(x) = —1, then we will have 
T n a(x) = +1 for all n > 1, hence <r is not an element of a stable pair, ii) Let C 6 {C e ,C°} and 
<T G 5o Suppose all the sites of A not belonging to the sea of chessboard are occupied by pluses, 
and suppose these pluses form a single cluster X C A. Consider the maximal Y C X such that for 
each y G Y there exists a 2 x 2 subset of Y containing y. If Y is not rectangular shaped, then there 
exists x such that a(x) = — 1 and at least two among its neighbors belong to Y. Then T n a(x) = +1 
for all n > 1 implies a is not an element of a stable pair. The proof can be easily generalized to 
the case with more than a cluster of pluses. Hi) The proof is similar to the one sketched for the 
case ii). □ 



Proof of Proposition 3.7, Let us consider a rectangle R^ m with 2<£<m<L — 2. 
Case 1: let 77 G <S_i be the trap such that VA\R e = ~-A\R e an< ^ r ^R l m = ^ > su PP ose ^ < A. 
Idea of the proof: we characterize the basin B(rj), that is we find T(r/); then we suppose the system 
prepared in 77 (namely o~o = v) anci > by means of Lemmata 13.41 and 13.51 we estimate both Tgyg^) 
and cr T „,-, , . 

By using (|3,19p . T(r/) can be estimated via the construction of the paths in H(r/). First of all 
we consider all the possible transitions that can be first steps for a path in 

1. the configuration on Ri yTn is flipped together with a minus spin in A \ Re >m adjacent to one 
of the plus spins in R^ m . A configuration r/i G S \ B(r/) is reached and #(77,771) — H(rj) = 
2(3(2 - h) =: $i. Hence 

T(tj) < H(n) + $1 = H (n) + 2/9(2 - /i) . (5.1) 

As shown in Fig. 13.31 the unique downhill path starting from 771 ends in the trap rj\ coincident 
with a chessboard inside a rectangle Rj/+i m (or 77J+1, depending on which side the protu- 
berance appeared) and with — 1 outside. We notice that the energy of 771 depends whether 
the plus protuberance appears in the middle of one side or on the corner, but H(r),r)i) does 
not depend on this detail. 

2. The configuration on i?£ >m is flipped together with a minus spin in A\R£ m with four neighbor- 
ing pluses. A configuration r)2 G B(rj) is reached such that H (77, 772) — H(rj) = 2/3(4 — h) > $1: 
this kind of steps can be neglected. 

3. All the spins inside Re m are flipped excepted one corner minus (if all the corner are pluses, 
then this step is considered after a full flip of the configuration inside Rt m ). A configuration 
773 G B(rj) is reached such that #(77,773) — H(n) = 2f3h. 

4. All the spins inside Ri^ m are flipped excepted one minus in the middle (not on the corner) of 
one of the four sides of the rectangle (if such a spin does not exist, this can happen in the 
case I = m = 3, then this step is considered after a full flip of the configuration inside Re )T n). 
A configuration 774 G B(rj) is reached such that #(77,774) — H(t]) = 2/3(2 + h) > <J>i: this kind 
of steps can be neglected. 

5. All the spins inside R# m are flipped excepted one minus with four nearest neighboring pluses 
(if such a spin does not exist, this can happen in the case £ = m = 3, then this step is 
considered after a full flip of the configuration inside Ri m ). A configuration 775 G B(rj) is 
reached such that #(77,775) — #(77) = 2/3(4 + h) > $1: this kind of steps can be neglected. 
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Fig. 5.6: Growth of a plus droplet inside the sea of minuses: appearing of a protuberance. 

6. All the spins inside R^ m are flipped excepted one plus spin. A configuration tjq G B{rj) is 
reached such that H{r],r]e) — H(rj) = 2/3(4 — h) > <&i: this kind of steps can be neglected. 

7. Two or more events among those listed above are performed simultaneously: the energy cost 
is smaller than 4>i only in the case of a simultaneous persistence of k minus corners of the 
chessboard. All the others multiple events can be neglected. 

From the list above it follows that there exists a path u/ = {rj, rji} G 3(7?) consisting of a single 
step of the first type; so ^ = H{rj) + <J>i = H(rj) + 2,3(2 - h). The only paths u G 3(r?) that 
can compete with u' are those whose first step is a single or a multiple minus corner persistence. 
After such a step (see, for instance, the first step in Fig. 13. 4p the configuration is a chessboard on a 
subset of R^m obtained by removing some of the four corners of i?^ m and —1 outside. By a direct 
inspection it follows that starting from this configuration the possible second steps of our paths are 
exactly those listed above: no new step enters into the game. 

By iterating the above argument, it follows that a path u" G 3 (77) such that $> w ii < $ w / 
can be obtained by using only single or multiple minus corner persistences. The best path lu" 
is a sequence of £ — 1 minus corner persistences on one of the two sides of the rectangle long £: 
<&w" = H{rf) + 2f3h{£ — 1). By comparing $> w ' an d arid recalling that £ < A, one obtains T(t/) = 
H(r]) + 2f3h{£ - 1). By Lemma S3 we obtain $(8(77)) = H(rj) + 2f3h{£ - 1) and U(B(r})) = W'}, 
where rj" is a configuration coincident with a chessboard in a rectangle i?^ m _i and with —1 outside. 
Finally, by applying Lemma 13.51 we can estimate "Tjy-g/ s ~ exp{2/3/i(£ — 1)} and we obtain that 
with high probability 0Y s ^g { } = Tj" . By using the Markov property and by iterating the argument 
above one completes the proof of part i) of Proposition 13.71 The proof of part ii) is similar. 
Case 2: Let us consider the trap r\ G S-\ such that Va\Ri = — -A\R e an< ^ = ^ ■ ^ ne 

proof is the same as in the Case 1. 

Case 3: Let us consider the trap r\ G 5_i such that il\\j^ e = — IaVR< anc ^ = ~^~—Tit ' ^-S am 
we suppose £ < A. As before we start by listing the transitions that can be first steps for a path in 
3(77): 

1. a minus spin adjacent to one of the four sides of the rectangle is flipped. A configuration 
r/i G S \ B(j]) is reached and H(r],r]i) — H{rj) = 2/3(2 — h) =: <&\. We notice that the unique 
downhill path starting from r/i ends in a trap rj as in Fig. 15.61 

2. A minus spin at distance greater or equal to two from any site of the rectangle is flipped. A 
configuration 772 G B(rj) is reached such that H(r], 7/2) — H{rj) = 2/3(4 — h) > $1: this kind of 
steps can be neglected. 

3. One of the four corners of the rectangle is flipped (a corner is a plus spin with two minuses 
among its nearest neighbors). A configuration 7/3 G B{rj) is reached such that ^(77,773) — 
H(ri) = 2/3/i. 
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4. One of the non-corner plus spin on one of the sides of the rectangle is flipped. A configuration 
774 £ B{rj) is reached such that -£/T(?7, 774) — H{rj) = 2/3(2 + h) > $1: this kind of steps can be 
neglected. 

5. One plus spin in the interior of the rectangle is flipped. A configuration 775 G B(eta) is reached 
such that H(r], 775) — H(rj) = 2/3(4 + h) > $1: this kind of steps can be neglected. 

6. Two or more spins are flipped simultaneously: the energy cost is smaller than $1 only in the 
case of a simultaneous flip of k corners. All the others many-spin flips can be neglected. 

From the list above it follows that there exists a path a/ = {77,771} £ ^(77) consisting of a single 
step of the first type; so $ w / = H(rj) + $1 = H{rj) + 2/3(2 - h). The only paths u G 3(77) that 
can compete with lj' are those whose first step is a single or a multiple corner erosion. Suppose 
that after the first step of our uphill path the configuration of the system is 773. Two more possible 
transitions must be taken into account in the analysis of the possible second steps: 

7. one corner spin is flipped. A configuration 777 £ B{~q) is reached such that H(rj^, 777) — if (773) = 
4/3/i. Indeed we have to take into account that the minus spin with two pluses among its 
nearest neighbors (the minus at the site flipped at the first step) must persist. 

8. The minus spin with two pluses among its nearest neighbors and one of its two adjacent plus 
spins are simultaneously flipped. A configuration 77s £ .6(77) is reached such that if (773, 77s) — 

H( m ) = 2/3/1. 

From the third step on no more possible transitions arise, excepted the obvious generalization of 
8: 

9. a corner plus spin at site x is flipped together with all the spins at sites y 7^ x such that 
Py(v y (u)^) — * 1 111 t ne limit (3 — > 00, where 77 denotes the actual configuration. The energy 
cost of this transition is 2/3h. 

We conclude that an estimate of T(??) smaller then = H (77) + 2/3(2 — h) can be obtained only by 
using an uphill path made of steps of types 3, 7, 8 and 9, or steps in which two or more transitions 
3, 7, 8, 9 are performed simultaneously. Consider a path obtained by using these transitions: until 
on each side of the rectangle there are two nearest neighboring pluses the configuration is still in 
B(i]). Hence, to exit B{rf) at least on one of the four sides of the rectangle there must be no pair 
of nearest neighboring plus spins. It is clear that the path u" made of steps 3, 7, 8 and 9, exiting 
/3(7/) and with minimal height along the path is the one described in Fig. 15,71 after a first step of 
type 3 and a second step of type 8, Z — 3 steps of type 9 are performed until the stable pair 77" is 
reached. The height along this path is <I> W » = H (77) + 2(5h(£ — 1). 

By comparing <I> W / and $> w ii and recalling that i < A, one obtains T(??) = = H (77) A-2j3h(l — 
1); by Lemma ETll we obtain ®(B(r])) = H(r]) + 2/3h(£- 1) and U(B{rf)) = {t]"}. Finally, by applying 
Lemma 13.51 we can estimate t s \^^ ~ exp{2f3h(£ — 1)} and we obtain that with high probability 
°' T s\B(r 1 ) = r l" • By using the Markov property and the results proven in the Case 1 one completes 
the proof of part i) of Proposition 13.71 The proof of part ii) is similar. □ 
The Propositions 13.81 and 13.101 can be proven via arguments similar to those used in the proof of 
Proposition 13.71 Proposition 13.91 is a byproduct of the proves of Propositions 13.71 and 13.81 
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Fig. 5.7: Shrinking of a droplet of pluses inside the sea of minuses. 
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